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Abstract. This paper is motivated by the current de- 
velopment of several space missions (e.g. ACES on In- 
ternational Space Station) that will use Earth-orbit laser 
cooled atomic clocks, providing a time-keeping accuracy 
of the order of 5 x 10~^^ in fractional frequency. We show 
that to such accuracy, the theory of frequency transfer 
between Earth and Space must be extended from the 
currently known relativistic order (which has been 
needed in previous space experiments such as GP-A) to 
the next relativistic correction of order 1/c^. We find that 
the frequency transfer includes the first and second-order 
Doppler contributions, the Einstein gravitational red-shift 
and, at the order a mixture of these effects. As for the 
time transfer, it contains the standard Shapiro time delay, 
and we present an expression also including the first and 
second-order Sagnac corrections. Higher-order relativistic 
corrections, at least 0{l/c^), are numerically negligible for 
time and frequency transfers in these experiments, being 
for instance of order 10~^° in fractional frequency. Par- 
ticular attention is paid to the problem of the frequency 
transfer in the two-way experimental configuration. In this 
case we find a simple theoretical expression which extends 
the previous formula (Vessot et al. 1980) to the next order 
1/c^. In the Appendix we present the detailed proofs of all 
the formulas which will be needed in such experiments. 



1. Introduction 

Recent advances in laser cooling of atoms have led to the 
development of a number of highly accurate atomic clocks 
(Caesium and Rubidium fountains) which have improved 
time-keeping accuracy by an order of magnitude during 
the nineties to currently sa 10^*''^ in fractional frequency 
(Laurent et al. 1999, Bize et al. 1999). For further im- 
provement of accuracy, the fountain clocks are limited by 
gravity, and therefore several experiments are planned for 
the near future that will fly laser cooled atomic clocks on 
board terrestrial satellites. One of these is ESA's Atomic 
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Clock Ensemble in Space (ACES) mission (Salomon & 
Veillet 1996) planned for 2005. This mission will place a 
laser cooled caesium clock together with a hydrogen maser 
on board the International Space Station. When com- 
bined with equally accurate time/frequency transfer sys- 
tems that allow the comparison between space and ground 
clocks, the experiment should become a useful tool for a 
number of applications in metrology, fundamental physics, 
atmospheric studies, geodesy etc. 

Therefore the ACES mission will include sufficiently 
stable optical and microwave time and frequency transfer 
systems to allow the comparison of the clocks with neg- 
ligible noise contribution from the transfer system itself. 
At the required accuracy this condition calls for two-way 
systems that exchange electromagnetic pulses in two direc- 
tions in order to eliminate or reduce a number of unwanted 
effects associated with the instrumental delays, the prop- 
agation within the ionosphere and troposphere, etc. The 
Time Transfer by Laser Light (T2L2) system (Fridelance, 
Samain & Veillet 1996) uses optical pulses that are emit- 
ted on the ground, reflected by the satellite and received 
back on the ground. The events of emission and reception 
are dated on the ground clock, and the event of reflection 
on the satellite clock. The microwave system is expected 
to exchange pulses in both directions and date all events of 
emission and reception. It will also measure the frequency 
of the emitted and received pulses on board and on the 
ground. Additionally it is expected to include two down- 
link signals at different frequencies in order to allow the 
cancellation of the residual ionospheric effect. 

At the required uncertainties a fundamentally rela- 
tivistic modelling of the experiment is indispensable. In- 
deed the relativistic effects when comparing two clocks 
separated by an altitude of « 400 km (as is the case for 
ACES) can amount to several parts in 10~*^, which ex- 
ceeds the expected uncertainties of the clocks (« 5x10^^'') 
by several orders of magnitude. Evidently, a relativistic 
treatment of distant clock comparisons using electromag- 
netic signals becomes necessary whenever the uncertain- 
ties of such experiments become smaller than the size of 
the relativistic corrections. This has already been the case 
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for the gravity probe A (GP-A) experiment (Vessot et al. 
1980), and is more generally the case since the advent of 
the Global Positioning System (GPS) which is now exten- 
sively used for distant clock comparisons. 

The first theoretical treatments of such comparisons in 
a relativistic framework were performed by Jaffe & Ves- 
sot (1976), Ashby & Allan (1979), Vessot, Levine et al. 
(1980), Allan & Ashby (1985), and Khoner (1993). Sub- 
sequent new methods with improved uncertainties, such 
as the Two Way Time Transfer (TWTT) (Hetzel & Sor- 
ing 1993) and the Laser Synchronization from Stationary 
Orbit (LASSO) (Veillet & Fridelance 1993), have led to 
the need for more accurate theoretical treatments that in- 
clude some higher order terms (Petit & Wolf 1994, 1997, 
Wolf k Petit 1995, Khoner & Fukushima 1994). The lat- 
ter papers treat in some detail the synchronization (time 
transfer) between distant clocks using electromagnetic sig- 
nals and the relation between the proper time of a clock 
on the Earth or on board terrestrial satellites and the geo- 
centric coordinate time, TCG. The relativistic theory for 
frequency transfer has been revisited recently by Ashby 
(1998), lea ding to results similar to those presented in 
Subsection |4.2| . In the perspective of the expected uncer- 
tainties of the ACES and other similar experiments, for 
time as well as frequency transfer, a re-examination of 
these formalisms has proved necessary. This is the subject 
of the present paper. 

Essentially we shall compute the relativistic transfers 
of time and frequency, including all the terms up to the 
order The coordinate time transfer up to this order 
is well known as it consists of the standard Shapiro time 
delay. Concerning the frequency transfer, we find that the 
formula for the one-way transfer is rather complicated. 
Our main result will concern the formula for the two- 
way frequency transfer, up to order 1/c^. This formula 
appeared previously in the paper by Ashby (1998) but 
without a detailed derivation. In this paper, we derive this 
formula and, more generally, we present a self-contained 
derivation of all the formulae needed in this context. These 
formulae will be of direct use in the ACES experiment, 
and a fortiori in more precise future experiments. Higher- 
order relativistic corrections are negligible with respect to 
the projected uncertainties associated with ACES. 

The numerical applications made in this paper con- 
cern the ACES mission with a transfer from the Space 
Station A orbiting at the altitude H = 400 km to a 
ground station located at B. For the velocities involved 
we use VA — 7.7 lO'^ m/s and vb — f ground — 465 
m/s; for the gravitational potentials, Ub/c^ — 6.9 10"^*^ 
and Ua/c^ = 6.510"^"; and for the Earth parameters 
GMe = 3.981014 mVs^ and Re = 6.37 10^ m. We con- 
sider that the experimental uncertainties of ACES will be 
at the level of 5 ps for time transfer and 5 x 10" for 
frequency transfer. 

The paper is organized in such a way that the main 
text summarizes all the results needed by an experimen- 



tal team in setting up a relativistic time and/or frequency 
transfer, with all the proofs and theoretical details rele- 
gated to the Appendix. In Section |2| we treat the problem 
of the transformation from proper time to coordinate time 
TCG at a level sufficient for ground clocks and the ACES 
space clock. In Section || we give the expressions required 
for time transfer (including Sagnac terms), and in Section 
^ for frequency transfer. Both the one-way and two-way 
transfers are considered in each case. 



2. Proper time in terms of coordinate time 

Throughout this work we use the geocentric inertial (non- 
rotating) coordinate frame GRS: Geocentric Reference 
System. Thus, x'^/c = t =TCG is the geocentric coor- 
dinate time[|, and x = (a;*) are the GRS harmonic spatial 
coordinates, for which the spatial metric is conformally 
fiat to order 1/c^. In these coordinates the metric interval 
including all the terms up to the order 1/c^ reads 



= _ ( 1 _ c^dt^ 



2U\ 

1 + ) S^jdx'dx^ , 



(1) 



where U is the total Newtonian potential, with the conven- 
tion that U >0 (lAU 1991), and where S.jdx'dx^ = rfx^ 
denotes the Euclidean space metric. The Solar-system 
barycentric coordinates, centred on the barycenter of the 
Solar system, are denoted T =TCB: solar system barycen- 
tric coordinate time, and X^: solar system spatial (har- 
monic) coordinates. The proper time of a clock A located 
at the GRS coordinate position TCA{t), and moving with 
the coordinate velocity v^i — d^iA/dt, is 



dTA 
dt 



= 1 - 



Li 
2 



C/e(xa) 

+V{Xa) - V{Xe) - x\d,V{XE) 
+x\Qi 



(2) 



Here, Ue denotes the Newtonian potential of the Earth at 
the position of the clock in the GRS frame, and V is 
the sum of the Newtonian potentials of the other bodies 
(mainly the Sun and the Moon), either at the position Xe 
of the Earth center of mass in barycentric coordinates, or 
at the clock location Xa- The three terms involving the 



1 The relation between t and the terrestrial time TT (re- 
alized by the International Atomic Time, TAI) is given by 
dt/dTT = 1 + Lg where Lg is a defining constant fixed in 
lAU resolution B1.9 (lAU 2000) as Lg = 6.969290134 10"^°; 
this value was chosen to be close to the previous definition 
Lg = Wo/c^ (lAU 1991) where Wo is the Earth potential (grav- 
itational plus centrifugal) at the reference geoid close to the 
mean surface of the oceans. 
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potential V represent the tidal field of the other bodies at 
the position of the clock. [To a good approximation (in- 
troducing errors of much less than 10^^^) one can write 
V{Xa) ~ V{Xe + xa)- Then, as usual, the tidal field at 
the position can be approximated using a Taylor ex- 
pansion by the standard expression ^x\x-'^dijV {X^) again 
neglecting terms smaller than 10~^^.] The terms due to the 
tides of the other bodies (i.e. involving V) are small for 
the ground station and the low orbit of ACES (of order 
2 X 10"^"^) and are either negligible or easily evaluated if 
required (depending on the final clock accuracies reached) . 
The last term in (^) is due to the non-geodesic acceleration 
of the center of mass of the Earth that is induced by the 
mass quadrupole of the Earth^ and which is numerically 
of the order of 3 x lO^^^m/s^. This term is negligible in the 
case of the envisioned experiments: numerically jQiX^/c^j 
amounts to less that 10"^". 

Thus, for application to ACES, and more generally 
to any experiment at a level of uncertainty greater than 
5 X 10^^^ on a satellite at similar altitude as ACES, we can 
keep only the first three terms in the relation (||) between 
the proper time ta and the coordinate time t: 



dTA 
dt 



Ue{-^a] 



(3) 



Note that at this level of uncertainty it is crucial to take 
into account in the relation (||) the non-sphericity (oblatic- 
ity) of the Earth Newtonian potential. In fact, it is even 
not sufficient to model the Earth potential with a J2-term 
taking into account the quadrupolar deformation. Rather, 
the potential U-e{'x.a) should be computed according to 
the detailed procedures of Wolf and Petit (1995), and Pe- 
tit and Wolf (1997). For example, for a clock B fixed on 
the Earth surface, the relativistic correction term appear- 
ing in the coordinate/proper time relation (^) is given with 
the required precision by 



2 pHb 

^ + C/e(xb) = Wo - / gdH , 
^ Jo 



(4) 



where Wq is the Earth potential at the reference geoid 
(Wo = 62636856 m^/s^), where g denotes the (gravita- 
tional plus centrifugal) Earth acceleration, and where Hb 
is the geometric height of the clock above the reference 
geoid. Note that this procedure is limited by the uncer- 
tainty in the determination of Wo which at present gives 
rise to an error of about 1 part in 10^^ in dTs/dt. For un- 
certainties below that level, it is expected that, inversely, 
clock comparisons with highly accurate space clocks will 
yield the best estimate of the potential on the ground. 

^ To a sufficient approximation we have 



where is the Earth quadrupole moment (see e.g. Brumberg 
& Kopejkin 1990 for details). 



3. Coordinate time transfer 

3.1. One-way signal transmission 

Let A be the emitting station, with GRS position Xyi(t), 
and B the receiving station, with position XB{t). We de- 
note by Ia the coordinate time at the instant of emission 
of a light signal, and by ts the coordinate time at the in- 
stant of reception. We put ~ |xA(iA)|, t^b = \xB{tB)\ 
and Rab — \^B{tB) — XA(tA)|, where | | is the Eucfidean 
norm associated with the metric 6ij . Up to the order 1 /c^ 
the coordinate time transfer Tab ^ tB ^ tA is given by 



R,4i3 , 2GA/E. fvA + rB+RAB 

Tab = h 5— In ■ — 

c c'^ \rA + rB- Rab 



(5) 



where the logarithmic term represents the Shapiro time 
delayQ (Shapiro 1964). See the Appendix for several 
derivations of the Shapiro time delay and for an alter- 
native expression given by ( A.40| ). In the case of zenithal 
geometry, i.e. propagation of the signal along the local ver- 
tical (for which jx^ — = |rB — r^l), between the orbit 
of ACES at 400 km (assumed in all numerical examples 
below) and the ground, the Shapiro time delay is 2 ps. In 
the case of zero elevation, it is 11 ps. 

In a real experiment, the position of the receptor B 
may be recorded at the time of emission tA rather than 
at the time of reception i^, i.e. we may have more di- 
rect access to XB(i^) rather than XB(tB), and the for- 
mula (^ gets modified by some Sagnac correction terms 
consistently to the order 1/c^. In this case the formula 
becomes 



Tai 



'Dab , D^B-VB(tA) 



c 

Dab 

2GA/E , 
^ — m 



(Dab.vb)' 

^ab 
rA + TB + Dai 
rA-\-rB -D 



Dab -as 



AB 



(6) 



where T)ab = ^B{tA) — x^(iA) is the "instantaneous" co- 
ordinate distance between A and B at the instant of emis- 
sion at A (we have D^^ — \^ab\), where VB{tA) denotes 
the coordinate velocity of the station B at that instant, 
and where as is the acceleration of B (in all the small 
correction terms of order 1/c'^ one can use, independent 
of the required order, the quantities at the instants tA or 
is). The second term in (^ represents the Sagnac term of 
order 1/c^ and can amount to 200 ns at low elevation; the 
third term or Sagnac of order is about 5 ps at low 
elevation (to be compared with the Shapiro term, which 
is 11 ps). 



In the case of a general metric theory of gravity, the factor 
two in front of the Shapiro time delay should be replaced by 
1+7 where 7 is the standard PPN parameter. 
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Fig. 1. Two-way time transfer in the non- rotating frame 



3.2. Two-way time transfer 

One signal is emitted from the satelhte A at instant tA 
and received by the ground station B at instant ts- A 
second signal (prime) is emitted by the ground B at in- 
stant tB' and received on A at tA' (see Figure]^). The two 
transmission times are 



Tab = tB ~tA , 

Tb'A' = tA' — tB' 



(7) 
(8) 



Furthermore, we denote the intervals of time between 
emission and reception on board the satellite A and at 
the station B by 



tAA' 
tB'B 



tA' — tA 
tB — tB' 



(9) 
(10) 



From these definitions we deduce the quantity At which 
is required for synchronization, namely 



At^tA- ti 



1 



{tB'B — tAA' + Tb 



Tai 



(11) 



Notice that A< is known because it is expressed in terms 
of the two transmission times Tab and Tb'A' which are 
computed from the theoretical formulas (|^) or (|^) valid 
for the one-way time transfer, and in terms of the time 
intervals tAA' and tB'B which are measuredon the satellite 
and on the ground respectivelj^. 

4. Frequency Transfer 

4-.1. One-way transfer 

The frequency transfer between two clocks requires the 
determination of the ratio f aI fB between the proper fre- 
quencies /a and /s delivered by the clocks on the satellite 

^ In the case of the time transfer by laser light T2L2 to be 
operated on ACES we have t^A' ~ since the signal is reflected 
instantaneously. 



(A) and on the ground {B). In practice this is achieved us- 
ing a transmission of photons from Ato B and the formula 



fA 
fB 



fA 
l^A 



VA 
VB 



VB 
fB 



(12) 



where va is the proper frequency of the photon as mea- 
sured on A (instant of emission tA)^ and vb the proper 
frequency of the same photon on _B at t s • The first bracket 
in (|l2|) is measured on A, the second bracket is given by 
the theoretical formula (13) below, and the third one is 
measured on B. In the one-way transfer of photons, as 
derived in the Appendix, we have 



1 1 

VA 


C/E(ri3)H 


2 




U^ivA) H 


2 



QB 



(13) 



For convenience in the notation, we henceforth denote the 
radial vectors (in GRS coordinates) by ta = XA{tA) and 
yb = x.B{tB)', and, as before, we have = {taI and = 
Irs I, as well as the coordinate velocities va — VA(tA) and 
vs = VB(iB). To the required order l/c^^, the last factor 
in (fl3) is obtained from 



= 1 



AGMe (rA + rs)NAij.VA 



Rab^-^ 



[rA + tbY 



r2 

^ab 



QB 



= 1 - 



AG Me (i^a + ?'s)Nab.vb - Rab' 



(r-A + r's)2-R; 



2 

AB 



(14) 



(15) 



with Kab rB - TA, Rab = |Rab| and Nab = 
Rab /Rab ■ See the Appendix for the derivation of these 
formulas. Note that the result (p^)-(p^) has been obtained 
in the Appendix assuming that the field of the Earth is 
spherically symmetric. Indeed, the J2-terms in the factor 
qa/qb do not exceed 4 x 10"^''. 

In the case of ACES, the various contributions in 
the one-way frequency transfer (p^)-(p^) are numerically 
as follows. First-order Doppler effect: for the satellite 
\Nab-^a/c\ < 2.6 X 10"^ for the ground |Nab.vs/c| < 
1.6 X 10^^. Second-order Doppler effect: v^/(2c^) 
3.4 X 10^1° for the satellite; v|/(2c2) < 1.3 x 10" 
for the ground. Gravitational red-shift (Einstein) effect: 
Ua/c^ = UE{rA)/c^ = 6.5 x 10"!°; Ub/c^ = 6.9 x 10-^°. 
The terms of order are less than 3.6 x lO^^"' for the 
satellite and 2.2 x 10^^''^ for the ground. 



4.2. Two- Way frequency transfer 

A "tracking" signal is sent from the ground station B at 
instant tB' , received on the satellite A at instant tA and 



< 
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used the upper bound (Rab)i 



tb/c^, where uj-e is 



the Earth's angular velocity of rotation. Now, according 
to (jl^), in order to reach the precision, the three 

previous terms are to be corrected by a factor that looks 
exactly like a first-order Doppler effect. Numerically, we 
have 



Fig. 2. Two-way frequency transfer in the non- rotating 
frame. 

instantaneously re-emitted by a satellite transponder to- 
ward B where it is received at instant ts- The down- link 
"clock" signal is emitted simultaneously with the track- 
ing signal at the transponding instant Ia, and received 
at tB (see Figure (Vessot et al. 1980). In the two-way 
frequency transfer the ratio va/vb in ([12|), needed for fre- 
quency comparison, is given by 



VA 2 \VB' 



^AB 



(16) 



The ratio {vb/vb') — {vBlvB')sta.t\o-a is measured at the 
ground station B, while I^ab is computed by means of the 
theoretical formula (derived in the Appendix) 



1 



MS 



Uab — 2^'ab ~ R-AB-as 



-v^.as + Hab-^b + 2vB.aB - wb-'^Ub) ■ 

(17) 



The difference of potentials between ground and satellite 
reads Uab = Ub — Ua where Ua = UE[TCA{tA)] and Ub = 
C/e[i"s(^s)] (note that Uab > 0); the gradient is VC/_b = 
dU-^i'T b) / dv B] the relative velocity is vab — v^(f^) — 
'VB{tB)', the acceleration of the ground is as — as (is); 
the derivative of acceleration is bs = da.B/dt. 

In the case of ACES we find the following numerical 
contributions. For the three dominant terms appearing at 
the order 1/c^, 



Ua 



Einstein 



^ 4.6 10-" 
1 o 



2d — Ol der Doppler 



2c2 



'AB 



A 



acceleration 



1 



T^AB-^B 



< 3.310-1" ; 
< 710-13 . 



(18) 
(19) 
(20) 



In the third term, or "acceleration" term, we have assumed 
that the station B is located at the equator, and we have 



< 2.710-^ 



(21) 



so this first-order Doppler factor induces a correction of 
the frequency shift at the minimal level 8.2 x 10-^''', which 
is measurable by ACES. Finally the four last terms in (17) 
are purely of order 1 /c^ . They amount respectively to the 
maximal values < 210-^^, 3.510-i^, and much less for 
the last two terms. These last four terms are in general 
negligible for ACES. The terms which are neglected in the 
formula (|l7|), which are at least C'(l/c'*), are numerically 
of the order of 10-^° or less, too small to be detected by 
ACES. 

In summary, the formula ( p^ and its derivation in the 
Appendix constitute the main results of this paper. To the 
dominant order l/c^ the expression was first derived by 
Vessot et al. (1980) who used it in their GP-A experiment. 
To the order l/c^ the formula (jl^) appears in a recent 
paper by Ashby (1998) but without detailed derivation. 
The final expression ( p^ is relatively simple, as compared, 
for instance, to the corresponding expression valid for the 

We find indeed 



one-way frequency transfer [see (13)-([15|) 
that many of the scalar products between ISSab and some 
velocities, though present in the intermediate steps of the 
calculation, drop out from the final result. The exception 
is for the factor in (|l^ which is made of the combination 
1 + ^ab-'^ab/c and which can be nicely interpreted as 
a modification, at the level 1/c'^, of the dominant term 
of order l/c^ by a first-order Doppler effect. The set of 
equations for time and frequency transfers given in this 
paper should be sufficient for the analysis of the planned 
clock experiments in Earth orbit, at the level 5 x 10-^^ in 
fractional frequency. 



Appendix A: Theory 

This Appendix presents several equivalent derivations of 
the formulas for the time and frequency transfers in the 
ACES experiment. Some of the basic material needed in 
these derivations is not new and can be found in standard 
textbooks such as Misner, Thorne & Wheeler (1973), and 
Will (1981). On the other hand, the problem of the propa- 
gation of light in a gravitational field has been solved in a 
general way at the linearized approximation: see Kopeikin 
& Schafer (1999) for a complete investigation and an entry 
to the literature. Here we use the explicit solution of the 
photon motion, the optical distance function for station- 
ary space-times (Buchdahl 1970, 1979), and the differen- 
tiation of the well-known Shapiro (1964) formula. 
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We shall first consider the transfer of coordinate geo- 
centric time t; for this purpose it is sufficient to approx- 
imate the gravitational field of the Earth as spherically 
symmetric (monopolar), and to neglect tidal terms, hence 
Ue = GM/r. However, as we have seen in Section||, higher 
spherical harmonics in the Earth potential are needed in 
the relation between coordinate and proper time. We de- 
note the mass of the Earth by M = Mb, the Carte- 
sian geocentric coordinates by {t,r), with t =TCG and 
r = X agreeing to this approximation with the geocentric 
spatial coordinates (r = |r|). Using spherical coordinates 
{t, r, 9, If} associated in the standard way to the geocen- 
tric coordinates (0<6'<7r, 0<(y9< 2tt), we can write 
the IPN metric as 

ds^ = ~f{r)c^dt^+g{r) [dr^ + r^d0^ + sin^ 0d(p^] ,(A.l) 
where to this order 



find immediately the first integrals kt = and = L, 
where E and L denote two constants along the trajectory 
(with i? > to ensure fc° > 0). Relating kt and k^ to the 
contravariant components of the wave vector, fc* = dt/dp 
and fc*^ = dip /dp, yields the two integrals of motion 



E = /(r)c- 



fir) 



2U 
72" 



1 



1 



2GM 
2GM 



(A.2) 
(A.3) 



As said above, the monopolar approximation to the Earth 
potential is sufficient for our purpose in most of this Ap- 
pendix. 

Photon equations of motion. In the geometric optics ap- 
proximation, where the photon's wavelength is much 
smaller than the typical size of the electromagnetic wave 
packet, as well as of the space-time radius of curva- 
ture, the photon's wave vector is null, fc^fc'' = 0, and 
is hypersurface-orthogonal, fc^ = dfj^S. The differential 
equation for the trajectory of the ray is k^ = dx^^/dp — 
g>"'{x)d^S{x). We have = Vf^ik^k") = 2fc'^V^fc^ = 
2k'^'Vi,k^, where we have used the fact that V^jfc^ is sym- 
metric as a consequence of fc^ = 9^5', so we find that 
k^^vkp, = 0, which is the geodesic equation, a more use- 
ful form of which reads 



dt 
dp ' 



L — g{r)r^ siv? 6 



dip 
dp 



(A.6) 
(A.7) 



dkf^ 
dp 



1 



k^ k'^ d f^g pa 



Here, p denotes an affine parameter along the trajectory, 
and the photon's null wave vector k^ is such that 



dx^' 
dp 







(A.5) 



(fc'^ is future directed, fc*^ — ck^ > 0). 

Although the metric (A.1)-(A.3) is spherically- 
symmetric, it is convenient to suppose that the motion of 
the photon takes place in an arbitrary plane, not necessar- 
ily the equatorial plane 9 = n/2. For instance, we can con- 
sider that 9 = represents the Earth's rotation axis, and 
that the emitted signal comes from a satellite moving on 
any orbit with inclination angle i with respect to the equa- 
tor. Since the metric is stationary and axi-symmetric, we 



Next, the equation corresponding to the 0-coordinate 
reads dkg/dp = g{r)r^ sin 9 cos 9 {k"^)^ . Inserting in this 
equation k'f = dp/dp as deduced from (A.7), and 
noting that g{r)r^ = ggg, we obtain d(kf)/dp = 
2L^ cos9 / sirC" 9d9 / dp. This equation shows that kg is a 
function of the ^-coordinate only, and we get k^ = r — 
L? I sin^ 9 where & is a new constant of the motion. Hence, 



(A. 



Finally, since g^^k^k^ = along a light ray, we get from 
the previous integrals of motion. 



9{r) E^ 
fir) c2 



62 



(A.9) 



It can be checked that the equation concerning the r- 
coordinate, i.e. dkr/dp = l/2kPk'^ d^gpa, is now auto- 
matically satisfied. Thus, the photon motion depends on 
the constants E, b and L; however, by eliminating the 
affine parameter p in favor of the coordinate time t we 
can parametrize the motion by only two constants, say 
b = bc/E and L = Lc/E. Finally, the solution reads 



(A.4) ^ 



dr 
cdt 

d9 

cdt 

2 d'P 
cdt 





9 




f r^ ' 




sin^ 9 


f L 




g sin^ 9 ' 



(A.IO) 

(A.ll) 
(A.12) 



We have introduced the signs e,- and eg of dr/dt and d9/dt. 
Because the metric (A.l) is spherically-symmetric, one 
could rotate the coordinates so that the plane of the mo- 
tion is simply the equatorial plane 9 = 7r/2. This would 
correspond to 6cquatoriai = L, and in that case pi would 
simply be the polar angle within the orbital plane. In our 
more general situation, it is not difficult to find the equa- 
tion of the orbital plane. We introduce the angle a defined 

by 



cos a 



&cos6' 

- L2 



sma 



/62 sin2 6l-L2 
62-L2 



(A.13) 
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In terms of a we can integrate the equation for the az- 
imuthal angle ip in the form 



tan((/3 — V'l) — "T tana 



(A.14) 



where (fi is an arbitrary constant. Clearly this is the 
equation of the orbital plane. Its line of node is de- 
fined by the direction cp — cpi — 7r/2 in the equatorial 
plane, and the inclination angle i of the orbit is such that 
sini = a/1 — L'^/b'^. The angle a is the polar angle in 
the orbital plane, oriented in the sense of the motion (we 
have da^ = dO^ + siv? Od(p'^). Introducing the Euclidean 
orthonormal basis e^, e^, e^p associated with the spheri- 
cal coordinates {r,9,ip}, the position and velocity of the 
photon read 



dx_ 

cdt 



^f{r) _ dr ^ 
g{r) cdt 



da 
cdt 



(A.15) 
(A.16) 



where Bq,, namely the unit vector in the direction of in- 
creasing a within the orbital plane, is given by 



£eee\ 1 



L2 



L 



62 sin^ e 



bsinO 



(A.17) 



In (A.16) we have introduced the unit tangent vector n 
along the photon's path, which is directly related to the 
contravariant components of the wave vector: = fc'/|k|, 
where |k| is the Euclidean norm, so that = SijU^n^ — 1. 
From the facts that fc'^ is a null vector and that the line 



element (A.l) is diagonal and spatially conformally flat, 
we can check that |k| = \/ f / 9 , from which we deduce 
that the covariant components ki are also proportional to 



n , VIZ 

ki 

ko 



9 k' 
lk° 



(A.18) 



From its definition ( A.16| ), and from the solution of the 
motion, we find that the tangent vector n is given by 



n = Sri 1 




(A.19) 



Next, the equations of motion for r{t) and a(t) in the 
orbital plane are 



dr 
cdt 
, da 
cdt 




(A.20) 
(A.21) 



The differential equation for the trajectory is 



dr 
da 



£rr\ 



/62 



(A.22) 



Photon trajectory. The trajectory at the relativistic order 
l/c2 (and even at order 1/c^), which is the solution of 
( A. 22] ), is an hyperbola whose focus is the center of the 
Earth, with impact parameter b and total deviation angle 
4:GM/c'^b. The equation of the path reads 



cos(q; — ao) = - 



2GM 
c^b 



(A.23) 



where the angle ao represents the direction of the peri- 
astron, at which the distance of closest approach is tq, 
which, to order l/c2, is given by 



ro = h- 



2GM 



(A.24) 



The tangent vector along the trajectory can be expressed 
in terms of the vectors e^g and erg , corresponding to the 
position (ro,ao) of the periastron: 

2GM 



c2ro 



sin(a — ao)er 



(A.25) 



As for the radial equation (A.20) it reads, after being ex- 
pressed in terms of r^ rather than b, 



cdt 



rdr 



1 



2GM 



1 



ro 



ro 



(A.26) 



This is readily integrated as 



Ct-tn = 



2GM 



r-ra 



In 



ro 



(A.27) 



where to denotes the instant of passage at the periastron. 
On the other hand, the trajectory is obtained as 



ao| — arctan 



ro 



2GM 
c^ra 



ro 



ro 



, (A.28) 



where ao is the angle at the periastron. 
Coordinate time transfer. We consider from now on a 
transfer from an emission point A along the trajectory 
to some reception point B. We introduce the Euclidean 
vectorial distance between the two points, defined as the 
difference of the coordinate positions of the points in the 
GRS coordinate system: 



'Rab — rs — ; R-As = |R.a_b| ; ^ab = R-ab/Rab ■ 

(A.29) 
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The position of the emitting point A is taken at the time 
of emission Ia, tca = ^^(iA), and similarly = YB{tB)- 
Notice that the norm Rab is simply the Euclidean norm; 
thus, Rab represents the "straight line" distance between 
A and B. In practice it is very useful to express all quan- 
tities in terms of such Euclidean notions of distance and 
direction as Rab and N^s- Once we have fixed our coor- 
dinate system to be (i, x), we can forget about the curved 



such that c{t — Ia) — V{TA^'r). This function is the so- 
called optical point characteristic of the stationary space- 
time, or optical distance between pairs of points (Buchdahl 
1970, 1979). In our case, we have 



V{va,v)^Ra 



2GM, frA 
— T\ — m 



Ra 



Ra 



(A.35) 



geometry ( A.l) and reason as if we were in flat space-time. 



Clearly, with the two points A and B being given, the 
trajectory is entirely determined. In particular, the radial 
coordinate tq of the periastron is uniquely fixed as a func- 
tion of r^, tb and Rab- Neglecting terms in 1/c^, we find 



where Ra = y — ya, Ra = Note that along the 

light ray from ya to r, the function S defined by S{t, r) = 
—ct + V{YA^Y) stays constant [S = —cIa)- Therefore S 
represents the phase of the signal and can be used to define 
the wave- vector a s — d^S. For this choice ko = —1 and 
ki = diV. Using ( A.35| ) we get 



4RiB 



AB 



[rA - tb) ]\{rA + tb) 



R 



AB 



AGM {rA + r)N\-RA'- 



[rA 



(A.36) 



(A. 30) where N^ = R^/R^- This is in complete agreement with 



The unit vector Nab is given as 



N 



AB 



2GM rs-TA 
( 

Rab 



(A.31) 



On the other hand, the unit tangent to the trajectory at 
the emission point A (say) is given from ( A. 25 ) as 



YiA 



2GM . , 

— s\a{aA - ao)ero 
c^ro 



(A.32) 



The difference between the two vectors ( A.31 ) and ( A.32 ) 
is a small quantity 0(l/c^). It is not diffi cult, w i th the help 
of our solution for the trajectory [cf ( A.23 )-( [A.25| ) and 
( A.3C )], to obtain the relation between these two vectors 
(always working consistently to the order 1/c^): 



Nab 
4GM 



Ras 



c^rA 



(rA + rBy-Ke 



YA 



iNAB-YA)^AB 

(A.33) 



The coordinate time transfer from Ato B, denoted Tab — 
tB — tA, follows from the equation (A. 27). By expressing 
it with the help of the previous notation, notably of the 
"straight line" distance Ras , we obtain the simple formula 



Tab — 



Rai 



2GM 



In 



TA 



TA + TB 



Rab/ 



(A.34) 



The second term, purely of order l/c'^, is the Shapiro log- 
arithmic time delay. 



The right side of (A.34) depends only on the spatial 
positions of the emission and reception points, i.e. ya and 
rs, and not on the instants Ia and/or ts- Indeed, for 
stationary space-times, the time t — tA elapsed from some 
emission instant Ia does not depend on tA but only on the 
spatial coordinates and r, so there exists a function V 



( [08| ) and ( IaJsI ) 



Given the simplicity of the result (A.34) for the time 



transfer when expressed in terms of the Euclidean distance 
Rab, one can guess that the formula can be derived di- 
rectly by integrating ds^ = along the path from A to 
B. Let us sh o w ho w this works. With the post-Newtonian 
metric (AJ_)-(A^) we have, along the photon's path, 



cdt= \ l 



2GM 



\dY\ 



(A.37) 



where jdrj is the Euclidean norm of the vector c?r = dx. 
Introducing Ra = y — ya as a new spatial coordinate 
along the path, we have \dY\ = {dR\ -|- R^dN^)^/^ where 
Ra = |Ra| and Na = Ra/Ra- But we know from ( A.33[ ) 
that Na differs from the unit tangent yia at the emission 
point by a small term 0(l/c^). So, when r varies (the 
origin point A on the path being fixed), we have ^Na = 
0{l/c^) and therefore we see that d'H\ — 0{\/c^) makes 
a negligible contribution. This shows that to order 1/c^ 
inclusively the time transfer can be calculated "along the 
straight line" ; we have 



dt= [l 



2GM\ dRA 
rc? I c 



and the total time transfer reads 



Tab — 



RAi 



2GM 



dRA 



(A.38) 



(A.39) 



To find the closed-form expression of the integral, we insert 
r = {R\ -t- 2RA.rA + which can be approximated, 

since Na — yia + 0{l/c^) and the integral already enters 
a small quantity, by r = (R^ -I- 2RAnA.rA + ^a)"^^^- Next 
we perform the integration over Ra from ^ to -B, and we 
are allowed to replace yia — Nab + 0{l/c^) within the 
result. Finally we find (see e.g. Will 1981) 



Tab — 



R 



AB 



2GM, /rs 
— ^ — In 



rs.N 



AB 



\rA + i"a-Nas 



(A.40) 
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This expression for the Shapiro time delay is shghtly dif- 
ferent from the previous form (A. 34) but can easily be 
reconciled with it. To this end one makes use of the iden- 
tities 



TA + ta-'Nab = 



jrB + Rab? -r\ 
2Rab 

^Rab 



(A.41) 
(A.42) 



which show that ( A. 34 ) and ( A. 40 ) are indeed totally 
equivalent. However, in practice, we shall prefer to use the 
formula ( A. 34 ) rather than ( A.4C| ) because of its structural 
simplicity. 



Sagnac terms. The formula (A. 34) gives the time transfer 
from the point ta = rA{tA) at the emission instant tA, to 
the reception point tb = YB{tB) at the reception instant 
Ib- But in fact, in the application to ACES, the useful 
time basis in the experiment is that provided by the clock 
at the point A (i.e., in the satellite) which records the 
instant of emission Ia- Therefore, it is more convenient to 



re-express the time transfer (A. 34) in terms of the position 
of the reception point B as it was at the instant Ia rather 
than at ts, i-e- when it was instantaneous with the instant 
of emission at A (in the GRS coordinate system). Thus, 
we introduce the "instantaneous" coordinate distance 



'Dab = YB{tA) - rA{tj 



D 



AB 



ID 



AB\ 



(A.43) 



and perform a consistent series expansion when 1 /c tends 
to zero. We know from ( A. 34 ) that the time transfer 
Tab = C(l/c)- Therefore, with the required accuracy, we 
can write Kab = D^s + TAB^B{tA) + ^TlgasitA) + 
0{l/c^), where is the acceleration. The successive rel- 
ativistic approximations are obtained by working out this 



formula iteratively together with (A.34). In this way we 
obtain 



Tai 



Bab , D^B-VB(iA) 



c 

Dab 
2GM 



■In 



{T>AB.^B? 

^AB 

rA + rB+ Bab 
ta + tb- Bab 



Das -as 



(A.44) 



where all quantities are measured at the emission instant 
Ia recorded by the on-board clock. The formula involves 
the Sagnac terms of first (1/c^) and second (l/c^) orders, 
as well as the Shapiro time delay [of course, consistently 
with the approximation, one can use indifferently in the 
Shapiro term either the distance Bab or RaBi and either 
TBitB) or rB{tA)]- 

One-way frequency transfer. In the geometric optics ap- 
proximation we have — d^^S and the phase differ- 
ence between two successive electromagnetic pulses is 



dS — {di^Sdx'^')A = {dfiSdx^)B, which in other words 
means {k^u^dT)A — {k^u^^dT)B where = dx'^/dr de- 
notes the four- velocity (such that g^^^u^u" = — 1). Now, 
dTA and drs are the proper periods of the same signal 
at A and B and we have va/vb = drB/dxA. Therefore, 
the one-way frequency shift of photons during the transfer 
from A to i? is 



i^A 

I'B 



{k^u'' 



{k^ut') 



(A.45) 



Separating out temporal from spatial components, we 



have fci, 



— ,,0 



u"(fco -I- kiv'^/c), where 



denotes the coordinate velocity and where u = 
{—gpav''v'^/c^)~^^'^- In a stationary space-time /cq is con- 
stant along the trajectory, i.e. (A:o)a = (A:o)b = —E/c in 
the notation ( [A.6| ); for instance, we have seen that with 
the choice S{t, r) = —ct -\- V{va^ r), where V is the optical 
distance ( |A.35| ), this constant is equal to minus one. Thus 
we can write 



VA 
fB 



(uO) 



l-f 



B 1 



(A.46) 



Next, we found in (A. 18) that the covariant components 
of the wave vector ki are proportional to the unit tangent 
vector n*. Using this we can further infer that 



VA 
VB 



/ON 1 _ , / gC'rA) ha.va 
{u°)a ^ V Kta) c 



(U°)i 



1 



firs) c 



(A.47) 



Within the first factor vP is the coordinate time vs proper 
time ratio, which has been computed in equation (H). 
Thus, from (m°)a/(u°)b = {dr / dt^)B / {dr / dt)A, we get 





C/E(rs) H 


2 




U^{va) H 


2 



(A.48) 



where Ua and Ub are the Newtonian potentials at the 
levels A and B. The factor (A.48) comprises the Einstein 



gravitational red-shift and the second-order Doppler ef- 
fects, which are both of order We recall from Section 
P that at the level of accuracy of ACES, the potentials in 
( [A.48| ) must take into account the oblaticity of the Earth. 



Now the second factor in the right-hand-side of (A.47) 
is nothing but the ratio of coordinate periods of the same 
signal at A and i?, namely dtB/dtA] it contains the first- 
order Doppler effect (~ 1/c) and the third-order (^ 1/c'^) 
terms we are looking for. For the computation of this fac- 
tor at the level 5 x 10~^^ we do not need to consider 
the J2 of the Earth potential, and we can approximate 
[/e = GM/r. In order to obtain a useful formula, we 
have only to substitute for the tangent vectors and 
wb their expressions in terms of the unit direction ^ab 
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into (A.47). The required relation between these vectors 
was found in (A. 33). Alternatively, one can use directly 
the formula (A.3(:). Finally, our end-result for the one-way 
frequency transfer with accuracy reads 



(A.49) 



where the ratio qa/qb is given by the formulas (p^-(p^ 
in the text. 

The previous expressions were obtained using our ex- 
plicit solution ( A.23)-( A.25) for the photon's trajectory. 
However, there is a simpler way to obtain the ratio qa/cib, 
which is based on the fact that since it is equal to the 
ratio of coordinate times: qa/^b = dtB/dtA, it can be 
computed directly by differentiating the coordinate time 
transfer Tab = ^b ~ ^a with respect to the emission time 
tA- The time transfer to order X/c" is given by the sim- 



ple formula (A.34) containing the Shapiro term. Thus, we 



must consider 



dTA^'--'^^ 

dtA 



2GM 



In 



TA + TB 



TA + rB - Rai 



. (A.50) 



The differentiation is to be performed taking into account 
the fact that the coordinate distance between A and B 
depends on both the emission and reception times, i.e. 
R.AS = \^B{tB) — ^A{tA)\- Thus we have for instance 



where (i^s/j^s')station is measured at the ground station, 
and where Aab is given theoretically by 



Aab = — 



VA 



1 

" 2 Vi/B 



(A.53) 



The ratios vb/va and vaIv'b are given by the formulas 
(|l3|)-(|l5|), valid for the one-way transfer, that we must 
simply insert into ( A.53 ) in order to get the result. The 
only problem is to correctly express all the quantities as- 
sociated with the up-link from the ground station B' at 
Ib' to the satellite AsXtA^ in terms of the same quantities 
associated with the down-link back from the satellite A to 
the station B at ts (see Figure |). Notably we need to 
express the unit vector N^/^ in terms of the unit vector 
Nab on the way back, as well as the velocity and accel- 
eration VB{tB) and asits) of the ground station at Ib- 
Neglecting 0(l/c^) we find 



Nb'A = - Nab< 1 + -Nab.vb 



[4(Nab.vs)' - 2v| - 2RAB.aB] 



[1 + -Nab-^tb) - ^RABaB ■ (A.54) 



We need also the velocity of the station at emission, 
v_B'(tB'); in terms of its velocity, acceleration and deriva- 
tive of acceleration at reception, VB(t_B), ^5(^5) and 
hBitB). We get 



dRAB f '^^B 

diA V diA 



(A. 51) V_B' = VB RABas + — [(R.AB-^B)aB 



We find that the ratio Qa/ib = dtB/dtA as obtained from 
the equation ( A.50|) agrees exactly with the one given by 
the expressions 

Two-way frequency transfer. In this method the "tracking" 
signal is emitted from the ground station at point B' and 
coordinate time tB', is received by the satellite at point 
A and time tA and immediately transponded back to the 
ground station where it arrives at point B and time tB- 
(A generalization of the formulae to the case where there 
is a time delay between reception and re-emission at A is 
also possible.) The "clock" signal is emitted at A, simul- 
taneously with the tracking signal when it is transponded 
at A, and transmitted to B. This technique, called the 
Doppler-cancelling technique, permits us to suppress the 
first-order Doppler effect and to drastically limit the un- 
certainties linked with the atmospheric contributions (see 
e.g. Kleppner, Vessot & Ramsey 1970). The ratio of the 
signal frequencies to be inserted in (|lj) follows from 



I^B 
VA 



VB 
VB' 



^AB 



1 

2 ' 



(A.52) 



R4 fib 



ab"b\ ■ 
(A.55) 



These relativistic expansions being fully taken into ac- 
count, we obtain finally, after a long computation, the 
formula (O) in the text, which is 



Aab = ^ 



2 n 

Uab — 2^-45 ^ R-AB-as 



Rab 



(-VA.as + RyiB.bB + 2v_B.a_B - vb-VUb) , 

(A.56) 



where Uab = Ub — Ua and v^s = — v^. AH the 
quantities at ^ or i? are expressed at the corresponding 
instants tA or respectively. 
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